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ABSTRACT 

We measure the mass of a modestly irradiated giant planet, KOI-94d. We wish to determine whether 
this planet, which is in a 22-day orbit and receives 2700 times as much incident flux as Jupiter, is as dense 
as Jupiter or rarefied like inflated hot Jupiters. KOI-94 also hosts at least 3 smaller transiting planets, all 
of which were detected by the Kepler Mission. With 26 radial velocities of KOI-94 from the W. M. Keck 
Observatory and a simultaneous fit to the Kepler light curve, we measure the mass of the giant planet and 
determine that it is not inflated. Support for the planetary interpretation of the other three candidates 
comes from gravitational interactions through transit timing variations, the statistical robustness of multi- 
planet systems against false positives, and several lines of evidence that no other star resides within the 
photometric aperture. We report the properties of KOI-94b (Mp= 10.5 ± 4.6Me, 7?p= 1.71 ± 0.16i?e, P = 
3.74 days), KOI-94c (Mp= IS.elfg^gM®, Rp= 4.32 ± 0.41i?e, P = 10-4 days), KOI-94d (Mp= 106± llMe, 
Rp= 11.27 ± 1.06i?®P = 22.3 days), and KOI-94e (Mp= 35^21^^®, Rv= 6.56 ± 0.62i?®, P = 54.3 days). 
The radial velocity analyses of KOI-94b and KOI-94e offer marginal (> 2a) mass detections, whereas the 
observations of KOI-94c offer only an upper limit to its mass. Using the KOI-94 system and other planets 
with published values for both mass and radius (138 exoplanets total, including 35 with Afp < 150M®), we 
establish two fundamental planes for exoplanets that relate their mass, incident flux, and radius from a few 

/ \ 0.53 / p \ -0.03 

Earth masses up to ten Jupiter masses: = 1.78 (^] I s-icm-2 ) for A/p < 150Af®, and 

/ \ -0.039 / \ 0.094 

Jg- = 2.45 [ff^] ( crgB-icm-2 ) > 150Af®. These equations can be used to predict the 

radius or mass of a planet. 

Subject headings: planetary systems — exoplanets, Kepler , stars: individual (KOI-94, KIC 6462863, 
2MASS J19491993-H4153280) — techniques: photometric — techniques: spectro- 
scopic 



1. INTRODUCTION 

One of the most pressing problems in exoplanetary 
physics is the anomalously large radii of close-in tran- 
siting gas giant planets. These hot Jupiters have radii 
larger than predicted by standard model s of giant planet 
cooling and contraction (for review s, see Fortney & Net- 
telmann||2010 Baraffe et al.||20l"0 ). Some mechanism, or 
a variety ot mechanisms, prevents planets from contract- 
ing, resulting in the observed inflated planetary radii. 
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The reasons for the radius anomaly for these planets 
could be tied to their formation and subsequent orbital 
evolution to close-in orbits, where the planets are subject 
to extremes of both tidal and radiative forcing. Mecha- 
nisms to explain the large radii of the planets have tried 
to tap the vast energy sources available in tidal or ra- 
diative forcing. It is critical to build up a large sample 
size of transiting giant planets, at a variety of orbital 
distances and incident fluxes, to better understand the 
p hysics that leads to the r adius anomaly. 



Miller & Fortney (2011) pointed out that all transit- 
ing giant planets receiving less incident flux than 2 x 10^ 
erg s~^ cm^^ do not appear inflated, meaning that they 
are all s maller in radius than ex pected for pure H/He 
objects. Miller fc Fortney] (2011 ) estimated the masses 
of heavy elements contained within t hat relatively cool 



sampl e of 16 planets (at that time). Demory & Seager 



(2011 1 extended this work to Kepler gas giant candidate 
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planets, and also found a lack of inflated candidates be- 
yond this same critical flux level, 2 x 10* erg cm"^. 
This incident flux is approximately equal to an equilib- 
rium temperature of 1000 K, for a zero Bond albedo and 
planet-wide redistribution of absorbed stellar flux. The 
detailed study of giant planets receiving less than 2 x 10* 
erg s"""^ cm~^ in incident flux will serve as a useful con- 
trast against the population of inflated giant planets that 
receive higher levels of incident flux. 

To probe the underlying physical processes that cause 
the observed diversity of planetary densities, we need to 
both expand our sample and to test links between the 
planets' physical properties and their orbital properties. 
The first part of the paper focuses on expanding our sam- 
ple. We measure the mass of a modestly irradiated giant 
planet or warm Jupiter, KOI-94d, in order to calculate 
its density and place constraints on its interior struc- 
ture. We wish to determine whether this planet, which 
is in a 22.3-day orbit and receives 2675 times as much 
incident flux as Jupiter (just a bit below the "critical" 
flux limit described above), is more similar to the bloated 
hot Jupiters or the cooler non-inflated gas giants, like our 
own Jupiter. 

In addition to the warm Jupiter, KOI-94 hosts at least 
3 smaller planets, all of which were detected through 
transit si gnatures in the phot ometry from the Kepler 
Mission. Hirano et al.l (2012) note the planet-planet 



eclipse that occurs in this system, which allows a de- 
tailed analysis of the planets' orbital dynamics. The 
multiplicity of this system presents an opportunity to 
examine the architecture of a closely packed system with 
a warm Jupiter. Using Keck-HIRES radial velocities 
(RVs), we measure the mass of the warm Jupiter. We 
obtain marginal mass measurements of two other plan- 
ets in the system and an upper limit to the mass of the 
fourth. Coupled with the mass of the giant planet we ob- 
tain from RVs, transit timing variations (TTVs) in the 
photometry allow an additional check of the RV masses 
obtained in this work. 

In the second part of this paper, we investigate how the 
mass a nd incident flux of a plan et relate to the p l anet's 



radius. Enoch et al. (2012) and Kane & Gelino (2012) 



have done similar work, but here we include many more 
low-mass planets (down to 3 M^), allowing us to probe 
the mass-radius-flux relation at lower masses than in ei- 
ther of those papers. 

This paper is structured as follows: §2 reviews mecha- 
nisms for inflating giant planets, §3 presents observations 
of KOT94, §4 describes the analyses used to derive planet 
masses, §5 argues for the planetary status of all four tran- 
siting candidates of KOT94 §6 describes the composition 
of KOT94d, §7 presents the radius-mass-incident flux re- 
lations for exoplanets and discusses possible interpreta- 
tions, and §8 summarizes the paper. 

2. INFLATION MECHANISMS 

A menagerie of radius inflation mechanisms have been 
proposed to explain the large radii of hot Jupiters. The 
most r ecent reviews, now be coming slightl y out of date 
are in |Fortney & Nettelmann (2010) and Baraffe et al 



( 2010 ). Without going mto detail here on any one mecha- 
nism, we classify possible explanations into three groups: 
incident flux-driven mechanisms, tidal mechanisms, and 
delayed contraction. 



Some inflation mechanisms are driven by incident flux 
from the parent star, also called insolation, whereby 
a small fraction of the absorbed stellar flux is trans- 
ported by a physical mechanism much deeper into the 
atmosphere, near or past the radiative-convective bound- 
ary. These mechanisms are discussed in a var iety of pa- 
pers: weather layer kinetic energy transpo rt (Showman 



,fc Guillot| 2002! 'Guillot fc Showman"2o62l), Ohmic dis 
sipation (|Batygin & Stevenson 2010; Pcrna et al. 2010), 
thermal tides ( Arras fc Socrates 2010), and mechanical 
greenhouse ( YoudrrTfc Mitclrcll 2010). These mechanism 
would in general aftect all close in giant planets to some 
degree, with the strength of the effect waning at lower 
insolation levels. Diversity in planetary mass, planetary 
heavy element masses, and planetary temperature would 
lead to a range of inflated radii. For reference, the time- 
averaged incident flux on a planet is: 



47ra2 
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(1) 



where i?* is the stellar radius, T^s is the effective stellar 
temperature, a is the semi-major axis, and e is the orbital 
eccentricity. 

Another class of solutions are tidal interactions be- 
tween the host star and planet, in particular, eccentricity 
damping. Coupled tidal and planetary structur e; evolu- 
tion has been calculated by a number of authors ( Boden- 



heimer et al.|2001 |Gu et al.|2003 


llbgui & Burrowsll^OOOl 


Miller et al. |2009| 


Ibgui et al.||20] 


L0| Leconte et al.|2010p. 



et al. (2010) who used the most detailed tidal evolution 
equations, is that since radius inflation via orbital ec- 
centricity damping is a transient phenomenon, it cannot 
be the "universal" radius inflation mechanism. Radius 
inflation by tidal heating is a short-lived phenomenon, 
but the average system age is several Gyr. However, in 
certain circumstances, including when a nonzero orbital 
eccentricity is maintained by outside forcing, tidal heat- 
ing can inflate giant planet radii for as long as the forcing 
lasts. The tidal power on a planet scales like: 



E : 



63 (27r)5 i?5 e2 

T G cKp^ 



(2) 



where E is the power, Rp is the planetary radius, G is 
Newton's gravitational constant, Q'p is the tidal dissipa- 
tion factor, an d P is the orbital period (rewritten from 
Equation 1 in Ibgui et al.| 2010 using Kepler's Third 
Law). Two advantages of rewriting this equation are 
that (1) orbital period, rather than semi-major axis, is 
an observable, and (2) this formulation of the equation 
does not depend on stellar mass. 

For completeness, we mention that delayed planetary 
contraction, due either to higher than anticipated atmo- 
spheric opacit ies ( Burrows et al.|2007|) or interior barrier 
to convection ( [Chabrier fc Barafte|2007 ) are another class 
of solutions which could contribute somewhat to larger 
radii. However, neither of these ideas should depend on 
proximity to the parent star, which should clearly be a 
feature of the correct solution(s). 

3. OBSERVATIONS 

In this section, we present photometry from the Kepler 
Space Telescope, as well as data from the ground-based 
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techniques of adaptive optics imaging, speckle imaging, 
and spectroscopy. Transits of the planet c andidates were 
identifi ed i n the Kepler l ight c urves in Batalha et al. 



(2012) and Borucki et al. 



( |2012| . In §4, we describe our 

method for simultaneously fitting the photometry and 
radial velocities derived from time-series spectroscopy. 
Here, we present adaptive optics imaging, speckle imag- 
ing, and spectroscopy of KOI-94. These observations rule 
out various instances of a nearby, stellar companion that 
could masquerade as a transiting planet, or possible false- 
positive scenarios, as discussed in §5. 

3.1. Kepler Space Telescope 

KOI-94 is on the edge of one of the chips of the Kepler 
CCD (see Figure [T]), so when the spacecraft rotates, the 
CCD loses the part of its field of view containing KOI-94. 
This causes two-quarter long gaps in the light curve. 




19.4 19.6 19.8 20.0 

RA (hours) 

Figure 1. A section of tiie Kepler footprint. Each point is a 
Kepler target. Tlie red five-pointed star marks tlic location of 
KOI-94. 

Nonetheless , the Kepler pipeline, as described in 
Batalha et al. (2012), has found 4 transiting planet can- 
didates associated with KOI-94 (see Figure [2]) . The light 
curve is phase-folded around their transit centers in Fig- 
ure ^ In summary, they are a super-Earth (KOI-94b) 
in a 0.7-day orbit, a mini-Neptune or super-Earth (KOI- 
94c) in a 10-day orbit, a Jupiter-size planet (KOI-94d) in 
a 22-day orbit, and a Neptune-size planet (KOI-94e) in a 
54-day orbit. The Kepler Input Catalog (KIC) ascribes 
an effective temperature of 6217 K and a radius of 1.238 
i?0to KOI-94, resulting in planetary radii of 1.41 J?®, 
3.44 i?0, 9.26 i?©, and 5.48 i?©for planets b, c, d, and e, 
respectively. Our analysis of the stellar spectrum, which 
finds different values for the stellar temperature and ra- 
dius and therefore for the planet radii, is described in 
§4. 

3.2. Adaptive Optics 

KOI-94 was observed with near-infrared adap tive op- 
tics on 2009 Nov 08 using ARIES on the MMT (lAdamsl 
et al. 2012). Images were obtained in both J and Ks, 
and reveal no companions closer than 7.5" (Aj — 2.5, 
Aifs = 2; see Figure [4|. The image FWHM was 0.23 in 
Ks and 0.43 in J. We can place a limit on undetected com- 
panions of A,7 = 2.2 and Aks = 3.4 at 0.5", Aj = 4.6 
and Aks = 5.9 at 1.0", and Aj = 8.7 and Aks = 9.1 
at 4.0" (and beyond). Any additional companions that 
would dilute the transit light curve and change the planet 



parameters are constrained to be faint or very close to 
the star. 

3.3. Speckle Imaging 

Speckle imaging of KOI-94 was obtained on the night 
of 19 June 2010 UT and the night of 23 Oct 2010 using 
the two-color DSSI speckle camera at the WIYN 3.5-m 
telescope on Kitt Peak. The speckle camera simultane- 
ously obtained 5000 (3000) 40 msec images on 19 June 
(23 Oct) in filters: V (center = 5620A, width = 400A) , R 
(center = 6920A, width = 400A) and / (center = 888OA, 
width = 400 A). These data were reduced and processed 
to produce a final reconstructed speckle image for each 
filter. Figure [5] shows the reconstructed R band image. 
North is up and East is to the left in the image and the 
"cross" pattern seen in the image is an artifact of the re- 
construction process. The details of the two-color speckle 
camera observations and the Kepler follow-up observing 
program are presented in Howell et al. (2011 1. 

On both occasions, the speckle data for this i?=:12.5 
star allow detection of a companion star within the ap- 
proximately 2.76 X 2.76" box centered on the target. The 
speckle observation can detect, or rule out, companions 
between 0.05" and 1.5" from KOI-94. The June 2010 
speckle image was obtained with the WIYN telescope 
during relatively poor native seeing near 1.0", while the 
Oct 2010 observations made during good seeing, 0.6". 
We found no companion star within the speckle image 
separation detection limits to a delta magnitude limit of 
^4 mag in the R band, 2 mag in the V band, and 3.4 
mag in the I band. 

3.4. Spectroscopy 

We obtained time-series spectroscopy of KOI-94 on the 
W. M. Keck I telescope with the HIRES echelle spec- 
trometer through an iodine cell. Most of the observa- 
tions occurred between May and August 2012, but our 
earliest observation was 2009-12-07. We rejected spec- 
tra with fewer than 4000 ADU (8760 photons), since 
these low-signal spectra resulted in large (> 6 m/s) ra- 
dial velocity errors. We also excluded spectra taken dur- 
ing a transit of KOI-94d to avoid confusion from the 
Rossiter-McLaughlin effect. This yielded 26 spectra that 
we deemed suitable for RV analysis, as presented in Ta- 
ble [1] We also obtained a spectrum of KOI-94 without 
the iodine cell as a template for RV analysis and for char- 
acterizing the star (see §4). 

Note that one observation, on JD 16135.095, occurred 
during a transit of planet e. To see whether this datum 
affected our mass results, we tried removing it and re- 
peating the circular analysis described in §4.2. Removing 
the datum changed the masses of the planets by ~ .Olcr, 
and so we can safely ignore the effects of this datum on 
our analysis of the planetary system. 

We determined the radial velocity of KOI-94 in each 
iodine spectrum based on the spectrum's Doppler shift. 
We used the lab- frame iodine lines that the iodine cell su- 
perimposes on the stellar spectrum to calibrate velocities 
to an instrumental precision of 3 ms~^ (although stellar 
jitter introduces additional errors) . This techniq ue is de- 
scribed in further detail in Howard et al. (2011). 



4. PLANET AND STELLAR PROPERTIES FROM RADIAL 
VELOCITIES, PHOTOMETRY AND SPECTRA 
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Figure 2. The top panel shows the aperture photometry long cadence (30 minute) light-curve. No corrections have been applied to the 
photometric measurements. The bottom panel shows the detrended light-curve. The data were detrended by applying a 2-day median 
filter. Observations that occurred during a planetary transit were excluded from the calculation of the median. The triangles indicate when 
transits occurred for planets b, c, d, and e with corresponding colours of cyan, green, red and blue. 

a circular solution. In our third and most sophisticated 
model, we fit for all four planets in eccentric orbits while 
simultaneously fitting the light curve. 



Table 1 

Spectroscopic Observations 



JD - 2440000.0 


Photons 


RV 


Error 








(m/s) 


(m/s; Jitter 


; 


15172.768 


9351 


-9.6 




6.1 


16076.070 


11709 


27.0 




5.5 


16100.047 


23095 


14.2 




5.0 


16109.929 


23065 


-13.5 




5.0 


16111.048 


22920 


-16.8 




5.0 


16113.043 


22594 


-14.9 




5.0 


16114.077 


14484 


-17.0 




6.0 


16115.074 


21584 


-5.0 




5.0 


16116.054 


23041 


8.1 




4.8 


16134.015 


22973 


-9.2 




4.9 


16135.095 


23592 


-9.7 




4.7 


16139.015 


22594 


6.8 




4.5 


16140.813 


22370 


3.8 




5.0 


16144.100 


20178 


20.3 




5.2 


16144.796 


23330 


14.6 




4.6 


16146.057 


18774 


18.5 




4.9 


16147.972 


22642 


2.5 




4.7 


16149.091 


15395 


-13.8 




5.0 


16149.752 


22918 


6.1 




4.5 


16149.766 


22629 


-6.7 




4.6 


16150.065 


19204 


-2.6 




4.8 


16150.079 


19162 


-5.6 




4.8 


16150.094 


18323 


-10.1 




4.9 


16151.075 


22791 


-10.0 




5.2 


16164.001 


23183 


23.0 




4.9 


16173.925 


21904 


-10.7 




5.1 



Note. — Observations with low signal-to-noise (fewer than 
8760 photons) were not used in our analysis and are omitted 
from this tabic. Photon counts arc averaged over all pixels. 

We use radial velocity measurements of KOI-94 to de- 
termine the mass of planet d, marginal mass detections of 
planets e and b, and an upper limit on the mass of planet 
c. A plot of radial velocity versus time and a four-planet 
circular fit is shown in Figure |6] 

We fit the RVs with three models, each of which pro- 
vides an interesting interpretation of the system. The 
first model has only one planet (the giant) in a circu- 
lar orbit. This is because the giant planet dominates the 
RVs, and so it is useful to compare a four-planet solution 
to a simpler one-planet solution. In our second model, all 
four planets are in circular orbits. Because KOI-94 is a 
closely-packed system, we do not expect large eccentric- 
ities of the planets, so we want to verify that a solution 
allowing eccentricities is not significantly different from 



4.1. Circular Orbit Solutions 

Here we compare the results of the one-planet and four- 
planet circular orbital solutions. For the one-planet fit, 
we find K = 16.3 ms~^, producing a reduced of 1.97. 
The data and phase-folded, one-planet circular fit are 
shown in Figure [7| 

In the four-planet circular model, the center-of-mass 
velocity 7 of the system and the semi-amplitude if„ of 
each planet n are allowed to vary, allowing 5 degrees of 
freedom; all other orbital parameters are fixed. We aug- 
mented the photon-noise errors by 3.0 ms~^in quadra- 
ture to account for the stellar jitter. The best four-planet 
fit had a reduced of 1.60. The best-fit radial velocity 
components from the four planets are shown in Figure [Sj 

Because the stellar jitter is unknown, we recalculated 
this fit, varying the stellar jitter to achieve = 1 (this 
yielded a jitter of 5.0 ms~^). The RV semi-amplitudes 
achieved in this fit were consistent with those assumed for 
a stellar jitter of 3.0 ms~^: the semi-amplitude for KOI- 
94d and KOI-94e changed by less than half a percent, 
and the semi-amplitudes for KOI-94b and KOI-94c fell 
within the la errors. We also note that because these 
are circular fits, the difference in between a model 
with a stellar jitter of 3.0m s~^ and 5.0m s~^ could arise 
from planetary eccentricities that are excluded from the 
model. Therefore, we assume a stellar jitter of 3.0 ms~^ 
and adopt the resulting parameters for the 4-planet cir- 
cular solution, which are reported in Table [2] 

We assessed the errors in Kn with a Markov Chain 
Monte Carlo (MCMC) analysis of 1 chain of lO'^ trials 
using the Metropolis-Hastings algorithm. Figure [9] shows 
the posterior likelihood distributions. The corresponding 
planet mass (Mp) and density (pp) distributions (which 
are calculated from the orbital period and planet radius) 
are also shown. These parameters are reported in Table 

m 

To calculate the probability of a non-detection in RVs 
for each planet, we assume a Gaussian posterior distribu- 
tion and calculate the fraction of trials that would have 
occurred for masses at or below zero, had we not im- 
posed a positive definite value for Kn- We calculate the 
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Figure 3. The top panel shows the photometric observations 
phased to the orbital period of planet b with the best fit transit 
model overlaid. The transit models for planets c, d and e have been 
removed. The transit times have been corrected for transit-timing 
variations. The next three panels show the transit lightcurves cen- 
tred on planets c, d, and e respectively. As in the upper panel, the 
best fit models for other transiting planet have been removed and 
corrected for transit timing variations. The colors of the overlaid 
models match the identification used in Figure [2] 



following probabilities of a non-detection via RVs: 0.04 
for planet b, 0.14 for planet c, 5 x 10~^^ for planet d, 
and 0.007 for planet e. However, these estimates of the 
non-detection probability are smaller than the true prob- 
ability of non-detection because of the Lutz-Kelker bias, 
especially for planet c. 

4.1.1. Limiting Outer Planets 



Table 2 




Planet parameters for circular orbits of KOI-94. 


Parameter 


Value 


Circular Keplerian Fit: KOI-94b 




1.76 it 1.4 


Orbital semi-amplitude K (ms^^) 


3.2 ± 1.7 


Mass Mp (M®) 


9.2 ±4.9 


Density pp (gcm~^) 


9.0 ±4.7 


Circular Keplerian Fit: KOI-94c 




Orbital semi-amplitude K (ms^^) 


1.6 ± 1.3 


Mass Mp (M®) 


6.5 ±6.3 


Density pp (gem"'') 


0.46 ± 0.37 


Circular Keplerian Fit: KOI-94d 




Orbital semi-amplitude K (ms^^) 


19.68 ±2.19 


Mass Mp (M^) 


102 ± 11.4 


Density pp (gem"'') 


0.380 ±0.042 


Circular Keplerian Fit: KOI-94e 




Orbital semi-amplitude K (ms~^) 


5.25 ±2.04 


Mass Mp (M^) 


36.6 ± 14.4 


Density pp (gcm~^) 


0.690 ±0.268 


Note. — The best fit parameters for four planets in 



circular Keplerian orbits, after adopting stellar parame- 
ters and orbital ephemerides from the eccentric solution 
described in §4.2 (see TablejsJ. 

The absence of a significant change in the radial ve- 
locity of the star between our earliest and most recent 
measurements strongly limits the possibility of a mas- 
sive outer companion. To quantify this, we computed a 
circular, four-planet orbit solution in which we allowed 
a linear trend in the radial velocity as a free parameter 
and used an MCMC analysis to explore the likelihood 
of a given linear velocity trend. Over a two-year base- 
line, the median trend was —0.0125 ± 0.0063 ms~^ per 
day, which corresponds to a 3g upper limit of -6.9 ms~^ 
per year. In the style of Winn et al. (2010), we com- 
pute the mass of an outer perturber based on the stellar 
acceleration, 7, assuming the planet induces Newtonian 
gravitational acceleration on the star and in the limit 
Mp <^ M*: 7 = GMp/a^. To induce a stellar accelera- 
tion 7 of -6.9 ms^^ per year via Newtonian gravity, an 
outer perturber would need to satisfy 



Mpsini 



(3) 



where i is the inclination of the planet's orbital plane 
with respect to the line of sight and Mj is the mass of 
Jupiter. Thus, with a significance of Su, we can rule out 
companions more massive than 3.9 Mj within 10 AU or 
more massive than 1.0 Afj within 5 AU. 

4.2. Eccentric Orbit Solution 

The four-planet fit in which we allow eccentricities to 
float is the most versatile model. This model has the 
advantage of simultaneously fitting the light curve and 
the RVs, which measures (thus refining A/^ and i?*). 
As demonstrated below, the values for planet masses in 
this model agree with the planet masses determined in 
the circular orbital solution to within Icr, and so we adopt 
the parameters from the eccentric solution for the rest of 
this work. 
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Figure 4. Adaptive optics image of KOI-94 (center) in J and Ks. The closest companion is 7.5" away. All of the objects in the 
frame are stars except for the spot in the bottom left, which is an artifact. 





Figure 5. Speckle image obtained 2012-10-23 at the WIYN tele- 
scope at 692 nm. The image spans 2.76" X 2.76". 
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Figure 6. Radial velocity versus time from May 2012 onward. 
Black points are data with la errors (assuming a stellar jitter of 
3.0 m/s); a circular four-planet fit is superimposed. 



Figure 7. A one-planet, circular fit to the RVs, phase-folded to 
the period of KOI-94d. The black points are the data (error bars 
are Ic), and the black line is the circular one-planet fit to the 
data. The gray points and fit are time-shifted repetitions of the 
black data points and fit. The red point is the oldest data point 
(2009); all other data are from summer 2012. 

The Kepler photometry and Keck radial velocities were 
simultaneously fit with an orbital model. The model has 
the following free parameters: mean stellar density (p*), 
scaled planetary radius (r„/i?*), impact parameter (6„), 
orbital period (fn), center of transit time (Tc^„), radial 
velocity amplitude (Kn), eccentricity (e„) and argument 
of pericenter via esiniUn and ecosuin. A photomet- 
ric and radial velocity zero point were also included. The 
number (n=l,2,3,4) corresponds to the parameters for 
planet b,c,d and e respec tively. The transit model uses 



the quadratic formulae of Mandel & Agol (2002). Limb- 
darkening coefficients were hxed in the m odels to 0.3236 
and 0.3052 as determined from the grid of Cla ret fc Bloe- 
(2011 1. The orbits are modeled with non-interacting 



men 

Keplerians. 

A best-fit model was initially computed by minimizing 
with a Levenberg-Marquardt style algorithm. This 
model was used to measure TTVs and to seed an MCMC 
analysis of the model parameter space. TTVs were de- 
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Figure 8. Radial velocity components from a four-planet circular fit. Each panel shows the radial velocity signature from one planet. The 
black line is the model fit; the black points are the RV data minus the model RVs from the other three planets. The blue points are the 
binned RVs. We chose the number of bins by rounding up the square root of the number of observations, creating 6 bins of equal spacing 
in phase. The error bars are the uncertainty in the mean of the data in each bin. The reduced of this fit is 1.60. The RVs provide a 9cr 
detection of KOI-94d, a 2.5<T detection of KOI-94e, a 2(T detection of KOI-94b and an upper limit on the mass of KOI-94c. 



termined by fitting for eacii individual transit, fixing all 
parameters except to their best-fit values. An up- 
dated best-fit model was then computed using the TTVs 
to produce a better phased transit for each planet. The 
time-series were corrected by computing time-corrections 
based on a linear-interpolation of the TTVs. 

Posterior distributions for each model parameter were 
determined with an MCMC-style algorithm . Th is model 
has been des cribed in|Gautier et al. ( 2012a[ ) and|Borucki 
ydifFcrcnce is that th 



hat the Tl'V mea- 



et al.| ( |2012[ ), the only 
surements are included and fixed to their best-fit values. 
Four Markov-chains were calculated each with a length of 
10^. The first 10% of each chain was discarded as burn- 
in. The median and ±la percentiles were calculated for 
each model parameter and reported in Table |3j which is 
in the appendix due to its length. 

4.2.1. Stellar Properties 

We used the template spectrum (without iodine) to 
determine the effective temperature (Teff= 6182 ± 58), 
surface gravity {\ogg= 4.181 ± 0.066), and metallicity 
([Fe/H]= -K0.6228± 0.0020) of KOI-94 through the Spec- 
tra Made Easy (SME) an alysis technique described in 
Valenti fc Piskunovj ( 1996| . Applying the evolutionary 
constraints of the Yi et al. (2001) model isochrones and 
the simultaneous solution of the lightcurve and RVs, we 
determine the mass {M^, = 1.277 ± 0.050) and radius 
{R^ = 1.52 ± 0.14) of KOI-94. The stellar properties are 
presented in Table |3j 



4.2.2. Properties of KOI-94d from the Eccentric Orbit 
Solution 

We detect KOI-94d with 9a confidence with the RVs. 
The eccentric orbit analysis gives a mass of 106 ± IIM^. 
This mas s is marginal l y con sistent with the mass re- 
ported by |Hirano et al.||2012[ ), who found Mp = 73± 25 
Me. 

4.3. Planet Masses from Transit Timing Variations. 

We observe coherent TTVs for KOI-94c and KOI-94d 
which are presented in Figure 11 TTVs usually indi- 
cate gravitational interactions between adjacent pairs of 
planets; such interactions allow us to refine the mass es- 
timates of these planets. 

We use the pres cription of Equations 8 and 9 from 
Lithwick et al.| f20"l2) to predict the anti-correlated TTV 
signals produced during the interaction of planets c and 
d between 650 and 800 days. Assuming zero eccentricity 
for both planets and ignoring factors of order unity, we 
calculate 



V = P 



7rj(2/3)(j - l)(i/3)A 



r = p' 



(4) 



(5) 



where ^ is the ratio of the inner planet mass to the stellar 
mass, /Li' is the ratio of the outer planet mass to the stellar 
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Figure 9. Posterior distributions for the four-planet circular fit. From left to right: likelihoods of K, Mp, and pp for (top to bottom) 
KOI-94b, KOI-94c, KOI-94d, and KOI-94e. The median, Itr, 2a, and Scr likelihoods are over-plotted in dark blue, cyan, green, and orange. 
To account for the positive definite bias of the likelihood distribution, a is measured above the median only; lower bounds are reflections 
of the upper bounds (i.e., a symmetric distribution is assumed). This enables a quantitative estimate of the number of trials consistent 
with a non-detection via RVs. Note that the best-fit values are obtained by minimizing x^; the posterior distributions are used only for 
estimating uncertainties. The best-fit values and uncertainties for the circular orbits are presented in Table [2] 



mass, V is the predicted semi-amplitude of the complex 
TTV signal for the inner planet, and V' is the predicted 
semi-amplitude of the complex TTV signal for the outer 
planet. A, which is given by 



A = 1 

P J 



(6) 



is the fractional departure from a, j '■ j — 1 mean motion 
resonance. 

Using the values for the orbital periods of planets c 
and d, for which the closest first-order mean motion 
resonance is 2:1, we calculate A = 0.074. Using the 
SME-determined mass of the star and the RV-determined 
planet masses, we calculate V = 7 minutes and V = 1.7 
minutes. These values agree with the observed TTV in- 
teraction in Figure [TT] for which the TTV amplitudes 
of planets c and d appear to be about 7 and 1 minutes, 
respectively. 

4.4. Dynamical Stability 

We investigate the long-term stability of the KOI-94 
system. We integrate the orbits of all four planets using 
the built-in hybrid symplectic/Bulirsch -Stoer integrator, 
part of the Mercury software package (Chambers 1999). 



We use the orbital and physical parameters of the four 
planets detailed in Table [3j The mean anomalies are de- 
rived from the best-fit joint radial-velocity solution. We 
define in the following a non-stable system when a close 



encounter between a given body and one of the four oth- 
ers occurs, within their common Hill radius. In a first set 
of 1-Myr integrations, fifteen eccentricity values for each 
planet are randomly drawn from the normal distribu- 
tion A^(e, (Tg), where e is the best-fit orbital eccentricity 
and CTe is the la eccentricity uncertainty. Eighty per- 
cent of these integrations yield close encounters between 
the two innermost planets b and c. Closer inspection 
of the results reveal that the proximity of b and c put 
stringent constraints on their orbital eccentricity, only 
allowing values less than or equal to the best-fit values 
for the system to be stable. 

In a second set of ten integrations, all planets are re- 
stricted to circular orbits with masses drawn from their 
normal distributions, in a similar manner as for orbital 
eccentricity in the previous step. All integrations with 
circular orbits resulted in stable systems, showing that 
varying planetary masses within their la uncertainty has 
a negligible influence on the dynamical stability of the 
system, if orbits are kept circular. 

More detailed dynamical analyses of the KOI-94 sys- 
tem would determine an upper limit on orbital eccen- 
tricity for the b and c components. This would however 
imply arbitrary assumptions on planetary masses, which 
would not completely exclude configurations with large 
orbital eccentricities. A possible resolution would be to 
precisely determine the planetary masses of the close-in 
components, either through additional RV observations 
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Figure 10. Radial velocity components from the four-planet ec- 
centric fit in which the photometry and radial velocities were fit 
simultaneously. Each panel shows the radial velocity signature 
from one planet (top to bottom: b, c, d, e). The red line is the 
model fit; the black points are the RV data minus the model RVs 
from the other three planets. The blue points are the binned RVs; 
their error bars the uncertainty in the mean. The shaded regions 
show phase-shifted repetitions of the data and fit. 

or an N-body analysis of the TTV signals from these 
planets. 

5. SUPPORT FOR THE EXISTENCE OF THE PLANETS 

In this work, we sought to measure the mass of the 
warm Jupiter KOI-94d to determine how the planet com- 
pared to the population of hot Jupiters. We simul- 
taneously fit the Kepler light curve and 26 RVs from 
Keck/HIRES to measure the mass of KOI-94d with a sta- 
tistical significance of 9a. We also sought to measure the 
masses of the other three transiting planet candidates in 
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Figure 11. Observed Transit Timing Variations in Q1-Q12 for 
(top to bottom) KOI-94b, KOI-94c. KOI-94d, and KOI-94e, with 
the same color scheme as Figures [2] and [s] The y-axis of each 
plot shows the difference between the observed transit time (O) 
and the transit time expected for a periodic orbit (C). Days are 
in JD— 2454900. The errors are smaller than the point size. We 
excluded one 5.5(t outlier in the O-C measurements for planet c at 
686 days. There is a section between days 650 and 800 in which the 
TTVs of planets c and d are anti-correlated, indicating a possible 
gravitational interaction. 

KOI-94 in order to better understand the architecture of 
this planetary system. Our RV measurements of KOI-94 
measure the masses of the other planets with signficances 
of less than 3a (see Tables [5] and [s]) . However, other ob- 
servations we have made support the interpretation of 
KOT94 as a system with four transiting planets. We 
describe support for the planetary interpretation of the 
candidates below. We defer our discussion of the com- 
parison of KOI-94d to other Jupiter-size planets to the 
next section. 

We report the properties of planets KOI-94b {Mp= 
10.5 ± 4.6Me, i?P= 1.71 ± 0.16i?e), KOI-94c (Afp= 
15.6l^5^6^e, Rp^ 4.32 ± 0.41i?e), KOI-94d (Mp= 
106± llMe, Rp= 11.27 ± 1.06i?e), and KOI-94e (Mp= 
astas^e- Rp= 6.56 ± 0.62i?e). Although the fidelity 
of the Ke pler candidates is very high (90 — 95% accord- 
ing to [Morton fc Johnson 2011 see Fressin et al. (in 



J If 

prep) tor an estimate ot hdelity as a function of plane- 
tary radius), false positives do exist among the Kepler 
planet candidates in the form of background eclipsing 
binaries, hierarchical triples, and other configurations of 
stars that result in the dilution of the eclipse signal, al- 
lowing it to masquerade as a planetary transit. The large 
size of Kepler pixels renders Kepler particularly vulnera- 
ble to bound companion stars with a large planet, such as 
a Neptune- or Jupiter-size planet, mimicking the transit 
signal of an Earth-size planet. To show that a plane- 
tary interpretation is superior to the interpretation of 
these various false positive signals, we outline the vari- 
ous measurements and statistical properties of the KOI- 
94 system that support the hypothesis of four transiting 
planets. 

5.1. Radial Velocities from KOI-94d 

The eccentric orbit solution of the RVs and light curve 
yield a semi-amplitude of 18.3 ± 1.5 ms^^ with a 22.3 
day period for KOT94d. Similarly, the four-planet circu- 
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lar orbit solution yields a semi-amplitude of 19.68 ± 2.19 
ms~^ to the RVs. The agreement between the circular 
and eccentric values for the semi-amplitude of K01-94d 
underscores the robustness of this measurement. 

At a 22-day period, a semi-amplitude of 18.3 ± 1.5 
ms^^ is consistent with the orbit of a planet around a 
star. A binary star system in 22-day period would have 
a velocity semi-amplitude of hundreds of kilometers per 
second. The width of the spectral lines is consistent with 
a stellar rotation speed of 8 kms~^, which is at least an 
order of magnitude smaller than the orbital speed of a 
such a binary system. Thus, we rule out the possibility 
of an eclipsing binary false positive in a 22-day orbital 
period for KOT94d. This supports the interpretation of 
KOI-94d as a planet. As discussed below, the planetary 
status of KOI-94d strengthens the case that the other 
transiting candidates are also planets. 

5.2. Observed TTV signature between KOI-94c and 
KOI-94d 

The apparent anti-correlation in the TTVs between 
planets c and d from 700-800 days (see Figure 11 sug- 



gests that these bodies are dynamically interactmg, i.e. 
that KOI-94c is also a planet. Our order-of-magnitude 
treatment of the TTV signatures indicates that the in- 
teracting bodies are of planetary masses, although the 
mass estimates are too low by a factor of 5. 

5.3. No Evidence of Another Star 

The adaptive optics images show no evidence of com- 
panions as close as 0.5" from KOl-94 within 2 magni- 
tudes in the J band or 3.4 magnitudes in the Ks band. 
Similarly, the speckle imaging shows no evidence of com- 
panions as close to KOl-94 as 0.6" within 4 magnitudes 
in the R band. Also, the 3a upper limit on an RV trend 
of 7 = 6.9 ms~^ per year between fall 2009 and fall 
2012 rules out a Jupiter-mass or more massive compan- 
ion within 5 AU. 

The detection of a second stellar spectrum in the spec- 
trum of KOI-94 would indicate background or compan- 
ion star. We searched for the spectrum of a second star 
in the iodine-free HIRES spectrum of KOI-94. To fit 
the spectrum of the primary star (KOl-94), we used a 
library of over 700 observed spectra from Keck HIRES 
that span T^g: 3266 K - 7258 K, \ogg : 1.46 - 5.0, and 
[Fe/H]:— 1.47 — hO.56. We found the spectrum from this 
library with the least squares difference from the spec- 
trum of KOI-94 (with similar results to the SME anal- 
ysis) and subtracted this spectrum. We then compared 
the residual spectrum to each spectrum in the stellar li- 
brary. The deepest minimum in between the residual 
spectrum and another star from the library was 2a at 
— 59.2kms^^; however, there were many other 2a so- 
lutions for that spectrum. We did not detect a sec- 
ond stellar spectrum with > 1% of the observed bright- 
ness of KOl-94 and a relative radial velocity of at least 
^ 8kms~^ with > 2a significance. This technique is 
sensitive only to neighboring stars within 0.4", the half- 
width of the slit of the HIRES spectrometer. 

5.4. Low False Positive Rate in Multi-Planet Systems 



Given the observed number of Kepler targets ut, the 
number of Kepler planet candidates ric, the number of 
Kepler multi-planet systems with i planet candidates 
nm.i, and the planet fidelity P (or single candidate false 
positive rate 1 — P), the fraction of systems with four 
planet candidates that we expect to consist of one false 
positive and three true planets (Pifp) is equivalent to 
the probability that a false positive is lined up behind a 
system that already has three true planets. The number 
of false positives among the candidate transiting plan- 
ets is (1 — P)nc, and the fraction of those aligned with 
apparent i-planet systems is Um^i/nt. Thus, 



{i-PW 



IFP 



{l-P)nc 
n-t 



(7) 



Adopting the values nt = 1 60171 from Lissauer et al. 
( [2012) and = 2300 from [Batalha et al. ( 2012 )7^Jid 
assum ing P = 0.9 (in accordance with Morton &: Johnson 
2011 1, we calculate 



IFP 



0.0014. 



The probability that 2, 3, etc. planet candidates are all 
false positives (un-associated eclipsing binaries that all 
happen to align behind KOl-94 within 0.5") is orders of 
magnitude smaller and can be ignored. The low false pos- 
itive probability is definitive of planethood even without 
the other arguments presented in this section. 

5.5. Rossiter-McLaughlin Measurement During the 
Transit of KOI-94d 



Hirano et al. (2012) measured the Rossiter-McLaughlin 
(RM) ettect durmg transit of KOT94d. They observed 
a clear RM signal that, when measured considering a 
transit depth of i?p/i?* = 0.06856 ± 0.00012, implied a 
projected stellar rotation of Vsmlg = 8.01 ± 0.73km s~^, 
which is in good agreement with their spectroscopically 
determined value of 7.33±0.32kms~^ and our SME anal- 
ysis. This co nstitutes eviden c e for the planetary status 
of KOI-94d. lAlbrecht et al.| ( |2013[ ) also measured and 
modeled the RM ettect during the same tr ansit of KOT 
94dand obtained results that agreed with [Hirano et al.| 
( [20T2| ). 

Modeling of the transit durations, including ingress 
and egress, of the KOI-94 transit candidates indicate 
small inclinations with res pect to the line of sight (see 



Lissauer et al. (2012) uses statistical arguments to cal- 
cuiate the false positive fraction in multi-planet systems. 



the inclinations in Table [3] Hirano et al. (j2012) measure 
a projected mutual inclination between planets d and e of 
S — —1.5° during a mutual transit event, indicating that 
these bodies are coplanar. That these planets are aligned 
with the stellar spin axis strengthens t he argument for 
their plan etary status. R ecent work by [Fabrycky fc Ke- 
[pTer Science Team (2012 1 shows that many Kepler multi- 
planetary systems are coplanar, so it is likely that planets 
b and c are coplanar with planets d and e. 

6. CONSTRAINTS ON THE COMPOSITION OF KOI-94d 

The goal of this work was to measure the mass of the 
giant planet, KOI-94d, and determine whether its bulk 
density was consistent with that of an inflated or a cold 
Jupiter. 

We have modeled the thermal evolution and contrac- 
tion of planet d using the methods described in |Fortney| 
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et al.| ( |2007[ ) and [Miller fc Fortney] ( |2011[ ). This model 
assumes no extra heat source irom the star. Including 
uncertainties in the system age, planet mass, planet ra- 
dius, orbital semi-major axis , and heavy element dis tri- 
bution within the planet (see Miller fc Fortney]|2011 ) we 
estimate that 18 ± 6 of heavy elements are contained 
within t he planet. This is very s imilar to estimates for 
Saturn ( [Saumon fc Guillot|[2004| . Based on the [Fe/H] 
of the parent star determined from our SME analysis of 
the spectrum, we estimate that the metals mass-fraction 
of planet d (^^pianct) is 11 ± 4 times that of the parent 
star. This metallicity enhancement, at this planet mass, 
ag rees well with other "wa rm Jupiter" planets studied 
by [Miller fc Fortney] ( f2011[ ). 

\Ve use state ot the art thermal evolution models for 
giant planets to establish that the bulk density of the 
planet is fully consistent with a non-inflated planet, and 
it is indeed enhanced in heavy elements compared to its 
parent star, in a manner similar to Saturn. We can fur- 
thermore compare KOI-94d to other planets by creat- 
ing mass-radius and mass-density plots for all planets 
with measured radius and mass (see Figure 12). These 
plots demonstrate the dependence of planetary radius 
(and density) on mass and incident flux. 

7. THE RADIUS-MASS-INCIDENT FLUX RELATION 

In this section, we examine empirical relations between 
the radius, mass and incident flux of exoplanets, includ- 
ing the KOI-94 system. We discuss possible physical in- 
terpretations of these relations and suggest avenues of 
future theoretical investigation. 

This pop ulation of planets w as compiled from exo- 
planets. org (Wright et al. 2011), which was queried on 
September 27, 2U12. Uur selection criteria were that the 
mass and radius of the planet were measured, and that 
errors in these measurements were reported, that the ef- 
fective temperature and radius of the star were measured 
with errors reported, and that the semi-major axis was 
measured with errors reported. This resulted in the sam- 
ple of 138 exoplanets listed in Table [4| 

The methods for determining the radius, mass, and 
incident flux of the planets in Table [4[ were as follows: 
All the planets transit, and so their raHius measurement 
was based on transit depth as determined within a self- 
consistent model of the observed light curve. For many 
of these planets, the uncertainty in planet radius is dom- 
inated by the uncertainty in stellar radius. The masses 
of the planets were measured by one of two methods: the 
majority were measured as Mpsini based on the radial 
velocity of the star; however, in s everal Kepler systems 



of multiple planets, TTVs aided (Cochr an et al. 2011 



Keple r- 18) or provided the sole means (Lissauer et al 



201 1[ Kepler-11) of planetary mass calculation. Oncer 
taint les in the planet mass stem from uncertainties in 
the stellar mass, uncertainties in the RV semi-amplitude 
K for low-mass planets, and uncertainties in the TTV 
analysis. The incident flux for each planet was calcu- 
lated using Equation [l] so uncertainties in the incident 
flux relate to uncertainties in stellar effective tempera- 
ture, stellar radius, and semi-major axis. 

There is a bre ak i n the mass-radius relation at ~150 
M(R (see Figure [T2[) . We chose 150 as the break 



based on a visual inspection of the mass-radius and mass- 
density diagrams. In determining the relation between 



planet radius, planet mass, and incident flux, we consider 
planets more or less massive than 150 separately. 
Our sample included 35 planets with Mp< 150 and 
103 planets with Mp> 150 M®. All four planets in the 
KOI-94 system are included in the low-mass population. 

To determine how incident flux affects radius, we cal- 
culated the time-averaged incident flux on each exo- 
planet from Equation [l] We divided the population 
into the "high incident nux" half (those with incident 
fluxes larger than the median incident flux, Fm = 8.6 x 
lO^ergs" cm~^), and the "low incident flux" half. These 



are shown in Figure 12 as the red (high flux) and blue 
(low flux) sets of points. For planets with Mp > 15OAf0, 
the planets that receive high incident flux are systemat- 
ically larger than planets that receive low incident flux. 

We performed a similar test to determine how the or- 
bital period affects radius (see Figure [l3[) . We divided 
the exoplanet population into those with lower than the 
median orbital period of 3.52 days (red, "short period") 
and higher than the median orbital period (blue, "long 
period"). Planet radius does not correlate with orbital 
period. 

Using the KOI-94 system and all other exoplanets with 
published values for both mass and radius, we establish 
two fundamental planes for exoplanets that relate their 
mass, incident flux, and radius from a few Earth masses 
up to ten Jupiter masses. We fit two planes between 
log(i?p), log(Afp), and log(i^), one in each mass regime. 
The resulting relations are: 



^ = 178 



0.53 



erg s ^ cm ^ 



-0.03 



(8) 



for Mp < 150M® 



and 



2.45 



Mp 



-0.039 



F 



erg s 1 cm ^ 



0.094 



(9) 



for Mp > 150Mff 



For completeness, we also flt two planes between 
log(/Op), log(Afp), and log(F): 



g cm 



Pp / Mp 



F 



0.09 



erg s 1 cm ^ ^ 
for Mp < 150Me 



(10) 



and 



g cm 



Pp ( Mp 

= 0.48 — ^ 



® 



1.10 



F 



-0.28 



erg s ^ cm ^ ^ 

for Mp > 150Me 



(11) 

Table [4[ lists the mass, radius, incident flux, orbital 
period, and reference for the planets used to calculate 
these fits. We include the Solar system planets in Figure 
[T2[for reference, although these planets were not used to 
generate the fits in Equations [8[-[TT| Slice through the 
planes at the median incident flux is shown as black lines 
in Figure [T2[ 

These planes were calculated by fltting a plane to the 
measurements of logMp, logi?p, and logF in each mass 



12 




100 
Mp/Mg 



10 7 



3 

°- 1 




100 1000 



Mp/M^ 



10 



cc 
tr 



1 Le 



H — + 




short period 
long period 



Figure 12. All planets with measured mass, radius, incident flux, and uncertainties therein, as listed in exoplancts.org fWright ct al. '20111 
see Table m. KOI-94 planets are plotted as five-pointed stars; Solar system planets are plotted as triangles. Left: Planet radius vs. planet 
mass. We divide the planets into two populations: those with higher-than-median incident flux (red), and those with lower-than-median 
incident flux (blue) The Solar system planets (purple) all receive less than the median incident fiux; KOT94 planets c, d, and e (cyan) 
receive less than the median incident flux, while KOI-94b (orange) receives more than the median incident flux. For Mp > ISOM^, higher 
incident flux correlates with larger planetary radius. For Mp < ISOAfg, higher incident flux {F) correlates with smaller planetary radius. 
The best-fit planes for Mp < 150JV% and for Mp > ISOM® are shown at the median fiux F = 8.6 X 10* ergs"l cm"^. Right: Planet 
density vs. planet mass. The coloring is the same as in the left panel, and the density fits for Afp < ISOM^and for Mp > ISOM^are shown 
at the median flux. For Mp > ISOMg), higher incident flux correlates with lower bulk density. For Mp < ISOAf^, higher incident flux 
correlates with higher bulk density. We determine empirical relations (see text) between log(Mp), log(i^), and each log(i?p) and log(pp) 
for Mp > ISOAfe and Mp < 150M^. 

[9] (i.e. for Mp > 150Me) were A = 0.39 ± 0.053, 
jB = -0.039 ±0.0096, and C = 0.094 ±0.0055. Thus, the 
dependence of radius on mass for low-mass planets is sig- 
nificant at lOcr, and the dependence of radius of flux for 
high-mass planets is significant at 17 a. The downward 
slope of radius versus mass for giant planets is detected 
at a significance of 4(t, and the downward slope of radius 
versus flux for small planets is uncertain. 

For Mp < 15OM0, the RMS scatter of the radius is 
1.41i?0 and the RMS scatter of the density is 2.69g cm~^. 
Considering that the average radius of a planet in this 
mass regime is 6.72i?0 and the average uncertainty in 
planet radius in this mass regime is O.34i?0(i.e. 5% of 
the typical planet radius), the RMS scatter of radii for 
Mp < I5OM0 is large compared to the uncertainties in 
measurements of planet radii. For Mp > I5OM0, the 
RMS scatter of the radius is 1.15i?0 and the RMS scat- 
ter of the density is 1.48gcm~^. Considering that the 
average uncertainty in planet radius for Mp > I5OM0 
is 0.76i?©, the RMS scatter is comparable to the uncer- 
tainty in planet radius in this mass regime. Interpreta- 
tions of the RMS scatter in each mass regime is discussed 
in §7.1. 

Note than for Mp < 15OAf0, radius depends strongly 
on mass (i?p cx Afp^^) and very weakly on incident flux 
(i?P cx F-°-"^). For Afp > 150M©, the dependence is 
reversed: Rp cx A/p*' '^^, and Rp cx i^O-C^. Since mass 
has little effect on radius for giant planets, the incident 
flux is the most important factor in predicting planet 
radius. 

In light of the very clear dependence of giant planet ra- 
dius on incident flux, and the possibility of a dependence 
of low-mass planet radius on incident flux, we wanted to 
examine the relations between incident flux a nd p lanet 
radius in greater detail. The top panel of Figure 14 shows 
planet radius as a function of incident flux for the low- 
mass and high-mass planets. The scatter in radius of 
the low-mass planets can be attributed to the strong de- 
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Figure 13. Planet radius vs. planet mass (same Figure \12\ but 
with a different coloring scheme). We divide the planets into 
those with lower-than-median orbital periods (red) and those with 
higher-than-median orbital periods (blue). The Solar System plan- 
ets (purple) and KOI-94 planets b, c, d, and e (cyan) all have longer 
than the median period. Orbital period does not correlate with 
planet radius. 

regime, with equal weight for each point. To see how 
the uncertainties in Mp, Rp, and F influenced the fit, 
we did 1000 trials in which we varied each measure- 
ment of Alp, Rp, and F based on a Gaussian distribu- 
tion with the la uncertainties reported for that planet 
to create a posterior distribution of coefficients. The 
median and average values for the posterior distribu- 
tion of the coefficients were consistent with the coeffi- 
cients of the original fit. If we write Equations [8] and |9] 
more generally as log(i?p/i?0) = A + Blog{Mp/A'I^) + 
Clog(F/erg s~^ cm~^), the la uncertainties in the coef- 
ficients for Equation Is] (i.e. for Mp < I5OM0) are A = 
0.25 ±0.185, B = 0.53 ±0.052, and C = -0.03 ±0.017. 
The la uncertainties in the coefficients for Equation 
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pendence of planetary radius on mass for Mp < ISOM^. 
However, the relation between the radii of giant plan- 
ets and the incident flux is clear in this plot because 
the dependence of planet radius on mass is very small 
for Mp > ISOMgj. 

It appears that for planets with Afy < 150M^, the 
planets receiving high incident flux are systematically 
smaller than planets receiving low incident flux. To ex- 
amine the validity of this correlation, we plotted the 
residuals to the relation Rp oc Mp-^^ as a function of 



incident flux (see the bottom panel of Figure 14). We 



found that the residuals only weakly depend on incident 
flux, but that there is a visible downward trend. Thus the 
suggestion that low-mass planets with high incident flux 
are smaller remains unclear. Future characterization of 
low-mass planets receiving varying amounts of incident 
flux will help elucidate this relation, if it exists. 
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Figure 14. Top: All planets with measured mass, radius, and 
incident flux, as listed in exoplanets.org (see Table 111. The black 
points are planets with Mp < 150Mg3, and the red points are 
planets with Mp> 150 Afg. For the giant planets, planet radius 
increases with incident flux. For low-mass planets, incident flux 
does not correlate with planet radius. This is because planet ra- 
dius scales with planet mass more strongly than with incident flux 
for low-mass planets, whereas incident flux is the primary factor 
in determining the radii of high-mass planets (see Equation [9]l . 
Bottom: Planet radius^ divided by dependence on planet mass 
according to Equation [s] versus incident flux for Mp < ISOA/g . 
The best fit to the data suggests a slight trend toward lower radius 
at higher incident flux, although the data are also consistent with 
no correlation between incident flux and planet radius. 



7.1. Interpretation of the Radius-Mass-Incident Flux 
Relations 

We consider the KOI-94 system in the context of other 
planets with published values for both mass and radius. 



We quantify these trends in two fundamental planes (see 
Equations [8]-|9| for exoplanets that relate their mass, 
incident flux, and radius from a few Earth masses up to 
ten Jupiter masses. These equations demonstrate that 
for low-mass planets, mass is much more important than 
flux in predicting a planet's radius, whereas for high-mass 
planets, incident flux is more important for predicting the 
radius of a giant planet than the planet mass. 

For low-mass planets, inverting Equation |8l predicts 
the planet's mass given its radius and incidentnux. The 
small coefhcient for mass for Mp > I5OM0 indicates that 
it is difficult to predict the mass of a giant planet given 
its radius and incident flux. 

The small RMS scatter in radius for giant planets 
{5Rp = 1.15 R^) compared to the typical uncertainty 
in radius (cr/jp = 0.76i?©) suggests that our model of 
mass and incident flux affecting the radius of giant plan- 
ets is appropriate and that other factors, such as orbital 
period and metallicity, play a small part, if any. 

For low-mass planets, the high RMS scatter {5Rp = 
1.41 i?©) compared to uncertainties in the radius (cr/jp = 
O.34i?0) indicates that additional physics might play a 
role in determining radius. In particular, the compo- 
sition of low-mass planets could strongly affect plane- 
tary radius. Planets on the low-mass branch could have 
compositions ranging fro m mostly hydrogen/helium (e.g. 
Saturn) to mostly water ('Charbonneau et al. 2009[ e.g. 
G J 1214 b) to mostly rock (Ijatal ha et al.|2Ull[ ""e:g^ Ke- 
pler 10b). The equations of state ^t these materials are 
quite different, allowing a planet of lOM© to vary by a 
factor of 5 in radius depending on its composition, in the- 
ory. Despite the potentially large range in compositions 
at a given mass, the low-mass fundamental plane works 
reasonably well: extrapolating this relation to Earth, we 
predict Rp = 1.15i?©. 

Planetary composition might be described in part by 
the mass-radius relation. For low-mass planets, the ex- 
ponential dependence of -Rp on Mp is higher than ex- 
pected. For a body of constant density (for instance, 
a rocky planet, if we ignore compressibility), we expect 
1/3 

Rp (X Mp . For low-mass planets, the observed rela- 



tion, Rp cx M^ 



.52 



is steeper. This steep increase in 
radius with mass cannot be explained by the compress- 
ibility of material, since compressibility would cause less 
increase in radius with increasing mass. The extra in- 
crease in radius per unit mass suggests a compositional 
gradient. Within the low-mass regime (Mp < ISOM©), 
higher-mass planets might have an increased admixture 
of volatiles. This is supported by observations in our own 
solar system; Uranus, Neptune and Saturn have a larger 
fraction of volatiles than Earth and Venus. 

In the giant planet population (Mp> 150 M©), the de- 
cline in radius with increasing mass corresponds to the 
onset of electron degeneracy as an important component 
of the planet's pressure along with Coulomb forces. For a 
body supported by electron degeneracy pressure, we ex- 
pect Rp (X Mp-i/3. The dependence of Rp cx Mp-° °'', 
which is measured with a significance of 4(t, indicates 
that Coulomb forces still play a significant role in sup- 
porting high-mass planets. 

In planetary modeling, 150 M© occ urs where the mass- 
radius relation for model planets (e.g. Fortney et al. 2007 1 
begins to gradually flatten out. It is marked by waning 
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relative importance of electrostatic forces which alone 
lead to i? ^ A/^/'^, and the gradual onset of degener- 
acy pressure, which for complete degeneracy leads to 
R ~ M~^/^. In an approximate way, 150 can be 
thought of as the start of the broad maximum in this 
curve that leads to radius being nearly in dependent of 
mass for gian t planets and brown dwarfs (jZapolsky &' 
Salpeter||1969[ ). 



The idea of searching for correlations between plane- 
tary radius, mass, incident flux, and other measurable 
planetary and stellar parameters is not novel; as men- 



i'he reversal of the correlation between Rp and F at 
150 is an interesting feature of the exoplanet popula- 
tion. The current population of observed planet s can be 
sculpted to some degree by evaporative mass loss ( Baraffe 
eF^ [20041 |Mi6l [Hubbard et aL||2007l [Lopez et a" 
2012p . X-ray and UV (XUV) photons can drive mass-loss 
ot hydrogen/helium planetary atmosph eres, and in th e 
energy-limited escape model (e.g., Erkaev et al. 2007), 
the mass loss rate depends inversely on the planet den- 
sity and linearly with the incident XUV flux. Planets 
with masses near ^150 are those that have low bulk 
densities (see Figure 12 ) are and generally more suscepti- 
ble to atmospheric mass loss. Under extreme XUV irra- 
diation, some rare planets may migrate from just above 
150 M® to below. 

Unlike the incident flux on a planet, the orbital pe- 
riod of a planet does not correlate with planet radius. 
This suggests that the mechanism that maintains the 
inflated radii of giant planets is driven by the incident 
flux rather than by tidal forces via eccentricity damping. 
However, because the eccentricities of many planets are 
poorly constrained, it is difficult to calculate the tidal 
power deposited in the planet. The uncertainty in the 
heat dissipation timescale further complicates our anal- 
ysis, since we cannot determine how long ago various 
planets might have been in sufficiently eccentric orbits 
for tidal forcing to inflate them. Regardless, a planet's 
radius is more strongly correlated with the incident flux it 
receives than its orbital period. Future studies of warm 
Jupiters and hot Jupiters with various orbital periods 
and eccentricities will help elucidate the role of potential 
interior heating mechanisms. However, the goodness of 
flt between the radius, mass, and incident flux for giant 
planets suggests that the current role of tidal heating 
(or any other inflation mechanism that is not driven by 
incident flux) is quite small. 

There are only 9 low-mass planets out of 35 that re- 
ceive more than the median incident flux and 9 that have 
shorter than the median orbital period, whereas the pop- 
ulatio n of gi ant planets is more evenly split. Howard 
et al. (2012) find that the occurrence of giant planets 
IS smaller tnan the occurrence of small planets at or- 
bital periods less than 10 days, suggesting that the hot 
Jupiters in this work are over-represented. The over- 
abundance of hot Jupiters compared to hot Neptunes in 
this sample could be due to historic observational bias 
of hot Jupiters, since the detection of Neptune-size and 
smaller planets from the ground was infrequent before 
the Kepler Mission. Measuring the masses of Neptune- 
size and smaller planets that receive high incident flux is 
necessary to probe the radius-mass-incident flux relation 
for low-mass planets. 



7.2. Comparison to Previous Work 



tioned in §1, Enoch et al. (2012) and Kane & Gelino 



( 2012 1 sought empirical relations between the properties 



ot planets and their host stars. Here, we incorporate an 
additional year's worth of planet discoveries and mass 
determinations, especially for Kepler planets. There are 
additional differences between our work and theirs, and 





2012 




m2 



The primary difference between Enoch et al. 
and the work presented here is that Enoch et al 
study 16 planets within orbital periods ot 10 days and 
with masses between 0.1 and 0.5 Mj, whereas here we 
study 35 planets with mass es below 150 M® (0.5 Mj). 
Whereas Enoch et al. (2012 1 consider three mass regimes 
of planetslD.l < Mp<jn>Mj, 0.5Mj < Mp < 2.0Mj, 
2.0Mj < Mp < 12Afj), here, we only consider two 
regimes (Mp < ISOM®, 150Me < Mp < 13Mj). With 
the additional transiting planets included here, we do not 
see evidence for three mass domains. 

However, we find that fits to the data incorporating 
only mass and incident flux p redict the radius of a planet 



as we l l as th e fits described in Enoch et al. ( 2012 ). Enoch 
et al. ( 2012 ) find an average absolute deviation m tne pre- 
dicted radiias from the true radius of 0.11 i?j (1.23 R^) 
across all three of their mass regimes, whereas we find 
a mean absolute difference of 1.41 for Mp < 15OM0and 
1.15 for Mp > 15OM0, or 1.23 i?®, over the whole sam- 
ple. Incorporating the orbital period of the planet, the 
stellar metallicity, and the stellar age does not signifi- 
cantly improve the accuracy of the predicted planetary 
radius. 

The slope of the mass-radius fit for Mp < ISOM® re- 
ported here, 0.53 ±0.005 2, is only 1.5g differe nt from the 
value of 0.45 reported in Enoch et al. (2012) for planets 
with Mp < 0.5Afj. Our slope ot the mass-radius fit for 
the high-mass planets, —0.039 , falls between the values 



for middle- and high-mass regimes in Enoch et al. ( 2012 ) 
and 



-0.09. 



Kane & Gelino (2012) used a similar prescription to 
our method to obtam a mass-radius relation. They fit 
a power law between the mass and radius of low-mass 
planets. However, they assumed that giant planet radius 
was constant with planet mass (at one Jupiter radius). 
Furthermore, they did not consider the effects of incident 
flux on planetary radius in either mass regime. 

8. CONCLUSIONS 

1. The KOI-94 system In this paper, we presented 
26 radial velocity measurements of KOI-94 ob- 
tained on Keck/HIRES. These measurements con- 
firm the giant planet KOI-94d and strongly support 
the existence of other transiting planets in the sys- 
tem. 

(a) Properties of Planet KOI-94d The mass is 
106±11M©. The radius is 11. 27±1.06i?©. The 
density is 0.363 ± O.lOlgcm"^. The planet is 
enriched in metals by a factor of 11 ± 4 with 
respect to the parent star. The mass of heavy 
elements, or "metals," in the planet is 18 ± 6 
Me. 

(b) Properties of Planets KOI-94b, KOI-94c, 
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KOI-94e These planets were detected at sig- 
nificance of less than 3cr in the radial velocity 
data. The radial velocity detections of planets e 
and b are marginal (> 2a), whereas the radial 
velocity measurements of planet c are consis- 
tent with a non-detection to la. More RVs and 
a numerical analysis of the TTVs are needed 
to better characterize these planets. However, 
the TTVs, multiplicity of the system, lack of 
evidence for another star, RM effect during the 
transit of KOI-94d and coplanarity of these ob- 
jects strongly suggest that these candidates are 
planets. 

(c) Dynamical Stability of the KOI-94 Sys- 
tem The system is dynamically stable on a 1 
Myr timescale for a variety of configurations, 
including circular orbital solutions. Although 
some eccentric solutions are stable, the best-fit 
solution considered in this work is unstable due 
to close encounters of planets b and c. 

(d) No Massive Outer Perturbers The non- 
detection of a trend in radial velocity, the lack 
of spectral features from a second star to 1% of 
the brightness of KOI-94within 0.4", and the 
non-detection of a companion outside 0.5"(in 
both AO and speckle imaging) rule out large, 
outer companions. 

2. Radius-Mass-Incident Flux Relation Using 
the KOI-94 system and other exoplanets (138 ex- 
oplanets total) with published values and uncer- 
tainties for planet mass, radius, and incident flux, 
we establish two fundamental planes for exoplanets 
that relate planet radius, planet mass, and incident 
flux between 2 and 3000 Af© in Equations [8]- 11 
The slope of the mass-radius relation for low-mass 
planets suggests that as low-mass planets increase 
in mass, the admixture of volatiles increases. Al- 
though the plane for for Mp > 15OM0 fits the giant 
planets very well, a higher RMS compared to uncer- 
tainties in radius ( 400%) for the low-mass planets 
suggests that additional physics, such as the com- 
position of heavy elements, might contribute to the 
radii of these planets. 
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APPENDIX 
LONG TABLES 

Table 3 

Star and planet parameters for the KOI-94 system. 



Parameter 



Value 



Notes 



Transit and orbital parameters: KOI-94b 



Orbital period P (days) 




3.743208 ± 0.000015 


A 


Midtrausit time E (B.JD) 




2454964.6175 ± 0.0021 


A 


Scaled semimajor axis a/R^, 




7.25 ±0.59 


A 


Scaled planet radius Rp/Ri, 




0.01031 ± 0.00014 


A 


Impact parameter b 




0.088 ± 0.072 


A 


Orbital inclination i (deg) 




89.30 ± 0.57 


A 


Orbital semi-amplitude K (ms' 




3.3 ±1.4 


B 


Orbital eccentricity e 




0.25 ±0.17 


B 


Center-of-mass velocity 7 (ms~ 




2.1 ± 1.4 


B 


Transit and orbital parameters: 


KOI-94C 




Orbital period P (days) 




10.423648 ± 0.000016 


A 


Midtrausit time E (BJD) 




2454971.00870 ± 0.00103 


A 


Scaled semimajor axis a/R* 




14.3 ± 1.2 


A 


Scaled planet radius Rp/Ri, 




0.02599 ± 0.00047 


A 


Impact parameter b 




0.41 ± 0.018 


A 


Orbital inclination i (deg) 




88.36 ± 0.75 


A 


Orbital semi-amplitude K (ms" 






B 


Orbital eccentricity e 




0.43 ±0.23 


B 


Transit and orbital parameters: 


KOI-94d 




Orbital period P (days) 




22.3429890 ± 0.0000067 


A 


Midtransit time E (BJD) 




2454965.74052 ± 0.00015 


A 


Scaled semimajor axis a/Ri, 




23.8 ±1.9 


A 


Scaled planet radius Rp/Ri, 




0.068016 ± 0.000080 


A 


Impact parameter 6 




0.055 ± 0.051 


A 


Orbital inclination i (deg) 




89.871 ±0.123 


A 


Orbital semi-amplitude K (ms~ 




18.3 ±1.5 


B 


Orbital eccentricity e 




0.022 ± 0.038 


B 


Transit and orbital parameters; 


KOI-94e 




Orbital period P (days) 




54.32031 ± 0.00012 


A 


Midtransit time E (BJD) 




2454994.2379 ± 0.0012 


A 


Scaled semimajor axis a/Rn, 




43.1 ±3.5 


A 
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Table 3 — Continued 



Parameter 

Scaled planet radius Rp / 
Impact parameter b 
Orbital inclination i (deg) 
Orbital semi-amplitude K (ms~^) 
Orbital eccentricity e 

Observed stellar parameters 

Effective temperature T(,ff (K) 
Spectroscopic gravity log(;(cgs) 
Metallicity [Fe/H] 
Projected rotation sin j(km s~^) 

Fundamental Stellar Properties 

Mass M*(Mq) 

Radius /{^(-Rq) 

Surface gravity log 3* (cgs) 

Luminosity L* (Lq) 

Kepler Magnitude Kp (mag) 

Age (Gyr) 

Planetary parameters: KOI-94b 

Mass Mp (M0) 

Radius Hp (iJ®) 

Density pp (gcm~^) 

Orbital semi-major axis a (AU) 

Incident Flux F (ergs"'^ cm~^) 

Equilibrium temperature Teq (K) 

Planetary parameters; KOI-94c 

Mass Mp (M®) 

Radius _Rp (i?®) 

Density pp (gcm~^) 

Orbital semi-major axis a (AU) 

Incident Flux F (ergs~^ cm~^) 

Equilibrium temperature Toq (K) 

Planetary parameters: KOI-94d 

Mass Mp (M®) 

Radius Rp (iJ®) 

Density pp (gcm~'^) 

Orbital semi-major axis a (AU) 

Incident Flux F (ergs"'^ cm~^) 

Equilibrium temperature Teq (K) 

Planetary parameters: KOI-94e 

Mass Mp (M®) 

Radius Hp (H®) 

Density pp (gcm~'^) 

Orbital semi-major axis a (AU) 

Incident Flux F (ergs"'^ cm~^) 

Equilibrium temperature Teq (K) 



Value Notes 



0.03960 ± 0.00024 A 

0.18 ±0.11 A 

89.76 ±0.15 A 

0.019 ±0.23 B 



6182 ±58 C 

4.181 ±0.066 C 

+0.0228 ±0.0020 C 

7.3 ±0.5 C 



1.277 ±0.050 D 

1.52 ±0.14 D 

4.181 ±0.066 D 

3.01 ± 0.60 D 

12.2 D 

3.16 ±0.39 D 



10.5 ±4.6 B,C,D 

1.71 ±0.16 A,B,C,D 

10.1 ±5.5 A,B,C,D 

0.05119 ±0.00067 E 

1.58 X lO'' A,C 

1486 F 



15.6I°5^g B,C,B 

4.32 ± 0.41 A,B,C,D 

0.9lt°'^? A,B,C,D 
0.1013 ±0.0013 E 
4.03 X 10^ A,C 
1012 F 



106 ±11 B,C,D 

11.27 ±1.06 A,B,C,D 

0.363 ±0.101 A,B,C,D 

0.1684 ±0.0022 E 

1.46 X 10* A,C 

806 F 



351^1 B,C,D 

6.56 ±0.62 A,B,C,D 

0.60l°'^^ A,B,C,D 

0.3046 ± 6.0040 E 

4.46 X 10* A,C 

584 F 



Note. — A: Based primarily on an analysis of the photometry, 

B: Based on a joint analysis of the photometry and radial velocities, 

C: Based on an analysis by D. Fischer o f the Keck/HIR ES template spectrum using SME ( [Valenti fc Piskunov|1996l ), 
D: Based on the Yale-Yonsei isochrones ( Yi et al. 12001 1 and the results from A, B, and C, 
E: Based on Newton's revised version of Kepler's i'hird Law and the results from D, 

F: Calculated assuming Bond albedos of 0.3 (b; Earth), 0.4 (c; Neptune), 0.34 (d; Jupiter), and 0.4 (e; Neptune) and complete 
redistribution of heat for reradiation. 



Table 4 

Exoplanets with Measured Mass and Radius 



Name 


Mp 


Hp 


Incident Flux 


Period 


First Ref. 


Orbit Ref. 




(Mffi) 


(fie) 


(ergs~^ cm^^) 


(days) 






55 Cnc e 
CoRoT-1 b 
CoRoT-10 b 


7.862 
327.284 
875.699 


2.078 
16.685 
10.862 


3.34E+09 
2.96E+09 
5.39E+07 


0.737 
1.509 
13.241 


iMcArthur ct al.'l'2004| 
Barge ct al. (::!U08f^~^ 
Bonomo ct al. {■ZUlO) 


Demory ct al.' 1*20111 

Barge et al. ( 200» ) 
IJionomo ct al. ('.iiUlO') 


CoRoT-11 b 


746.274 


16.013 


2.04E+09 


2.994 


'Gandom et al. (:^OiU) 


Uandom ct al. (20iU) 
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Table 4 — Continued 



Name 




Rp 


Incident Flux 


Period 


First Ref. 


Orbit Ref. 




(Me) 


(Ke) 


(ergs~^ cm~^) 


(days) 
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Table 4 — Continued 
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